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,  ABSTRACT 

p- 

An  efficient  contention  resolution  scheme  has  been  developed  by 
Gallager  to  be  used  in  allocating  usage  of  a  common  channel  to  a 
large  number  of  Independent  transmitters.  Time  is  divided  into  equal 
size  intervals  called  slots-  Message  transmissions  begin  only  at 
slot  boundaries  and  do  not  overlap  boundaries.  Depending  on  the 
channel  history,  as  monitored  by  the  users,  the  algorithm  defines  a 
system-wide  ^transmission  interval*  at  the  beginning  of  each  slot. 

A  user  can  only  transmit  when  it  possesses  a  message  with  a  generation 
time  falling  in  the  current  transmission  interval.  This  algorithm 
is  modelled  as  a  Markov  process. 

The  Gallager  algorithm  is  based  on  the  assumption  that  the 
transmitters  make  no  errors  in  the  detection  of  channel  activity. 

This  thesis  investigates  the  case  where  a  limited  class  of  detection 
errors  is  introduced  into  the  system.  A  modified  algorithm  is 
proposed  to  effectively  correct  these  errors.  This  Noisy  Channel 
Algorithm  includes  a  stack  mechanism  and  two  special  processing 
states.  The  algorithm  is  again  modelled  as  a  Markov  process  and 
retains  some  key  features  of  the  Gallager  algorithm.  The  algorithm 
is  analyzed  and  the  degree  to  which  system  performance  is  downgraded 
by  errors  is  determined. 
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CHAPTER  I 


I.  INTRODUCTION 


1.1  Description  of  the  Problem 


This  thesis  is  concerned  with  the  problem  of  efficiently  allocating 
the  use  of  a  noisy  channel  to  a  large  number  of  independent  transmitters. 
The  transmitters  are  attempting  to  communicate  with  a  common  destination. 
Simultaneous  transmission  by  two  or  more  channel  users  is  termed  a 
"conflict",  and  results  in  the  central  receiver  being  unable  to  detect 
any  coherent  message.  A  retransmission  scheme  must  therefore  be  estab¬ 
lished  that  will  resolve  conflicts  and  provide  high  system  throughput. 

A  first-come  first-served  (FCFS)  algorithm  to  accomplish  this,  based 
on  the  assumption  of  error-free  detection  when  no  conflict  occurs, 
has  been  developed  by  Gallager.  This  algorithm  will  be  presented  and 
then  used  as  the  basis  for  designing  a  more  general  algorithm  that 
allows  for  a  limited  class  of  detection  errors. 

In  this  first  chapter  the  source  and  channel  models  are  specified, 
and  a  few  common  accessing  strategies  are  briefly  outlined. 


1. 2  Source  and  Channel  Model 

In  modelling  this  communication  system  we  first  assume  that  there 
are  an  infinite  number  of  transmitters,  or  message  sources.  Each 
message  is  of  a  fixed  length,  called  a  packet.  The  time  necessary  to 
transmit  a  packet  is  defined  as  a  slot.  All  time  quantities  referred 
to  in  this  paper  will  be  in  terms  of  this  fundamental  unit,  the  slot. 


The  users  generate  messages  in  a  Poisson  manner  with  a  global  rate 
of  X  packets/slot.  Thus  the  probability  of  a  particular  source  gene¬ 
rating  a  message  is  0,  but  the  expected  number  of  messages  generated 
system-wide  during  each  slot  is  X.  At  most  one  message  can  be  success¬ 
fully  transmitted  during  each  slot,  and  it  is  immediately  apparent 
that  X  must  be  less  than  one  to  insure  system  stability. 

The  users  are  not  capable  of  communicating  with  each  other,  but 
they  can  listen  to  the  channel  and  instantaneously  sense  whether  0,  1 
or  more  than  one  (a  conflict)  messages  are  being  transmitted.  If  any 
portion  of  two  messages  overlap  during  transmission,  both  messages  must 
be  retransmitted  in  their  entirety.  We  begin  by  assuming  that  receiver 
and  user  detection  is  error-free.  Thus  if  only  one  source  is  using 
the  channel  its  message  is  received  correctly  at  the  central  facility. 
Different  classes  of  system  errors  will  be  explored  later. 

It  is  important  to  note  that  messages  are  generated  continuously 

by  the  sources,  but  are  only  sent  when  the  source  determines,  based 

on  the  algorithm  implemented,  that  it  should  transmit.  This 

algorithmic  decision  is  a  function  of  the  generation  time  of  the 

message  and  the  history  of  channel  transmissions  as  detected  by 

the  source.  In  evaluating  a  conflict  resolution  scheme  we  want  to 

* 

determine  the  greatest  arrival  rate  X  such  that  the  expected  number 

* 

of  messages  generated  but  not  yet  sent  remains  bounded.  This  A  is 
called  the  capacity. 

1.3  Some  Accessing  Strategies 

Two  commonly  employed  multi-accessing  strategies  are  time  division 


multi-accessing  (TDMA)  and  frequency  division  multi-accessing  (FDMA) . 
Because  we  are  concerned  with  a  large  number  of  users,  FDMA  is  expensive 
in  terms  of  hardware  costs.  The  bursty  nature  of  the  sources  being 
considered  leads  to  poor  delay  characteristics  with  a  TDMA  approach. 
Although  they  prevent  the  occurance  of  conflicts,  neither  TDMA  nor 
FDMA  is  particularly  well  matched  to  efficiently  allocating  a  channel 
to  a  large  number  of  bursty  users. 

Another  class  of  algorithms  that  have  been  analyzed  is  the  ALOHA 
system,  developed  at  the  University  of  Hawaii.  In  this  strategy,  when 
a  conflict  occurs  the  contending  sources  wait  a  random  length  of  time 
before  retransmitting.  A  variation  on  this  theme,  that  increases 
capacity  from  ^  to  — -  , is  called  slotted  Aloha.  Slotted  Aloha  places 
the  constraint  that  all  transmissions  must  begin  at  discrete  times 
separated  by  one  slot  intervals.  Both  of  the  Aloha  schemes,  however, 
become  unstable  as  the  number  of  users  goes  to  infinity. 

Several  strategies  have  been  proposed  that  involve  using  a  secondary 
channel  to  reserve  time  on  a  primary  communication  channel.  The  secondary 
channel  is  randomly  accessed  and  the  primary  channel  is  allocated 
through  dynamic  TDMA  to  those  sources  requesting  it.  This  concept  is 
based  on  the  assumption  that  any  reservation  message  is  much  shorter 


than  a  packet.  A  system  of  this  type  again  runs  into  stability  problems 
as  the  number  of  sources  gets  very  large.  Also,  this  approach  still 
requires  an  Aloha  -  type  algorithm  (to  access  the  secondary  channel) , 
and  is  really  not  offering  any  new  contribution  to  the  fundamental 
multi-accessing  question. 
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Another  Aloha  variation,  carrier  sense  multiple  access,  requires 
users  to  listen  to  the  channel  before  transmitting  to  ascertain  if 
there  is  a  transmission  already  in  progress.  This  is  only  useful  when 
the  time  necessary  to  listen  to  the  channel  is  small  in  comparison  to 
a  slot.  With  many  users  this  algorithm  has  stability  problems  induced 
by  many  sources  trying  to  transmit  upon  the  termination  of  a  previous 
message . 

In  a  1977  Ph.D.  thesis  [1],  Capetanakis  proposed  a  fundamentally 
different  approach  to  resolving  conflicts.  He  begins  by  assigning  an 
address,  corresponding  to  a  terminal  node  on  a  binary  tree,  to  all 
sources.  When  a  conflict  occurs  during  a  slot,  only  the  sources  in 
one  half  of  the  tree  are  permitted  to  retransmit  during  the  next  slot. 
Successive  binary  division  of  the  sources  continues  until  the  conflict 
is  resolved.  The  algorithm  then  works  back  toward  the  root  node,  level 
by  level,  by  transmitting  all  sources  in  the  branches  that  had  been 
inhibited.  If  a  conflict  needs  to  be  resolved,  successive  binary 
divisions  again  take  place.  When  the  root  node  is  reached  all  conflicts 
from  the  original  transmission  have  been  resolved.  It  is  important  to 
note  that  all  messages  generated  during  the  process  of  resolving  a 
conflict  cannot  be  transmitted  until  the  return  to  the  root  node.  Thus 
no  messages  are  sent  for  the  first  time  during  the  conflict  resolution 
phase.  Capetanakis  later  generalized  his  algorithm  by  allowing  the 
tree  root  node  to  be  of  degree  greater  than  two,  depending  on  the 
length  of  the  previous  conflict  resolution  period,  and  thereby  achieved 
a  capacity  of  .430. 


The  Capetanakis  algorithm  is  important  in  that  it  geometrically 


decreases  Che  number  of  sources  able  to  transmit  during  successive 
slots  of  conflict  resolution.  A  binary  elimination  takes  place.  This 
is  as  opposed  to  TDMA  or  FDMA  schemes  that  linearly  allocate  the 
channel  or  schemes  where  all  sources  randomly  access  the  channel 
(Aloha) . 

This  concept  of  systematically  inhibiting  the  sources  allowed  to 
transmit  underlies  the  Gallager  algorithm.  While  Capetanakis  classifies 
sources  by  a  static  binary  address,  however,  Gallager  looks  at  only  the 
sources  that  have  generated  messages  and  classifies  them  in  terms  of 
the  time  their  message  was  generated. 

In  the  next  chapter  the  Gallager  algorithm  and  its  significant 
features  are  explained,  as  well  as  some  proposed  modifications  to  the 
algorithm.  In  addition  a  Markov  process  model  for  the  Gallager  algorithm 
will  be  presented.  The  Markov  model  is  essential  to  the  analysis  of 
the  performance  of  this  algorithm  and  a  modified  model  will  eventually 
be  used  to  represent  the  Noisy  Channel  Algorithm. 

II.  A  CONFLICT  RESOLUTION  ALGORITHM 
2. 1  The  Gallager  Algorithm 

The  Gallager  algorithm  is,  like  the  Capetanakis  tree  algorithm, 
based  on  a  slotted  system.  By  this  we  mean  that  transmissions  can  only 
begin  at  discrete  times  (beginning  at  a  slot).  We  assume  that  all 
channel  users  are  synchronized  with  the  slot  boundaries.  The  key  idea 
in  the  Gallager  procedure  is  that  at  the  beginning  of  each  slot  a  system- 
wide  transmission  interval  is  defined.  During  that  slot,  only  sources 
that  have  messages  with  generation  times  falling  in  the  transmission 
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interval  are  permitted  to  transmit.  If  a  conflict  occurs,  the  next 
transmission  interval  will  be  the  first  half  of  the  previous  transmission 
interval. 

In  the  Gallager  algorithm  the  messages  being  processed  during  a 
given  slot  were  actually  generated  at  an  earlier  point  in  time.  Thus 
at  each  slot  there  is  a  system  lag.  If  the  system  is  to  be  stable  this 
lag,  on  the  average,  must  not  be  increasing  with  time. 

In  implementing  the  Gallager  algorithm  each  source  must  keep 
track  of  three  variables.  The  system  lag  -  d,  the  length  of  the  trans¬ 
mission  interval  -  8,  and  the  generation  time  of  any  message  that 
particular  source  has  not  yet  successfully  transmitted.  At  the  start 
of  a  slot,  all  messages  generated  in  the  interval  begun  d  slots 
ago,  will  be  sent.  The  transmission  interval  length  and  lag  are  updated 
system-wide  after  each  slot  on  the  basis  of  that  slot's  result  (0,1  or 
more  than  one  message  sent) . 

The  workings  of  this  algorithm  are  best  illustrated  by  going 
through  an  example.  In  figure  1.1  the  system  is  in  a  renewal  state. 

By  this  we  mean  that  all  messages  generated  before  the  point  A  have 
been  successfully  transmitted,  and  nothing  is  known  about  the  message 
distribution  beyond  A  (except  that  it  is  Poisson).  The  policy  at  a 
renewal  state  is  to  transmit  an  interval  of  length  8.  beginning  at 
the  point  A.  8q  is  to  be  chosen  so  as  to  optimize  system  throughput. 

If  d  <  8q  at  a  renewal  point,  we  assume  that  the  system  either  trans¬ 
mits  an  interval  of  length  d,  or  waits  until  d  >  8  before  transmitting. 

—  o 

Note  that  in  real  time  we  are  at  the  slot  beginning  at  B  when  this 
first  transmission  interval,  beginning  at  A,  is  sent.  In  figure  1.1 
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there  is  only  one  message  transmitted,  and  thus  the  system  progresses 
to  another  renewal  state.  We  note  that  one  slot  has  elapsed.  The  system 
variables  are  updated  as  follows: 


d 

new 


£ 

new 


dold 

l 


+ 


1  - 


l 

o 


If  the  channel  had  been  idle  during  this  slot  the  algorithm  would  have 
continued  identically.  Again,  all  messages  in  the  interval  (zero  in  this 
case)  would  have  been  successfully  transmitted. 

In  the  next  slot  transmission  (fig.  1.2a)  a  conflict  occurs. 

We  define  this  interval  as  a  conflict  interval.  The  algorithm  states 
that  in  resolving  a  conflict  we  define  the  next  transmission  interval 
as  the  first  half  of  the  conflict  interval.  Hence  in  1.2b, 


d  *  1  +  d  .  , 
new  old 

l  .  .  * 

0  _  old  _ 


new  2  2 

During  this  slot  we  find  that  the  channel  is  idle  and  we  therefore 
know  that  there  are  >  2  messages  in  the  second  half  of  the  conflict 
interval.  The  first  half  of  the  conflict  interval  has  been  processed 
and  we  define  a  new  conflict  interval  equivalent  to  the  second  half  of 
the  original  conflict  interval.  Again  we  want  to  send  the  first  half 
of  this  conflict  interval  (the  third  quarter  of  the  original  interval) 

So  in  1.2c, 

d  ■  d  ..  +  1  -  £  .  . 
new  old  old 


old 


new 


o 
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The  next  transmission  results  in  a  conflict.  At  this  point  an 


important  feature  of  the  Gallager  algorithm  comes  into  play.  If  the 
first  half  of  a  conflict  interval  contains  a  conflict,  the  second  half 
of  the  interval  is  returned  to  the  waiting  interval.  The  waiting  Interval 
is  defined  as  the  set  of  times  that  will  not  be  part  of  a  transmission 
Interval  before  the  system  has  passed  through  at  least  one  renewal 
state.  In  figure  1.2d  the  waiting  interval  begins  at  C  and  continues 
to  the  present  time. 

In  returning  a  portion  of  the  conflict  interval  to  the  waiting 
interval  we  are  making  use  of  a  property  of  the  Poisson  distribution. 

For  a  Poisson  distribution  of  messages  over  an  interval,  the  nunfter  of 
messages  in  each  of  two  disjoint  portions  of  the  interval  are  independent 
random  variables.  And  for  two  independent  random  variables,  x^  and  x2» 

Pr(x2  -  x|xj  2,  xx  +  x2  2)  *  Px  (x) 

In  other  words  the  interval  being  returned  has  a  Poisson  distribution 
of  messages.  This  preserves  the  fact  that  any  portion  of  the  waiting 
interval  has  a  Poisson  message  distribution. 

At  the  next  slot  (1.2d), 


l..  f. 

5  .  .  ° 

new  2  8 

There  is  a  new  conflict  interval  and  its  first  half  is  being 
transmitted.  A  successful  transmission  results  and  we  are  left  with  an 
interval  that  contains  1  or  more  messages  (1.2e).  The  algorithm  states 
that  in  this  instance  we  send  the  entire  interval,  regardless  of  its 


length.  Hence 


No  conflict  occurs  during  this  slot  and  we  have  reached  a 
renewal  state.  The  original  transmission  interval  of  length  has 
been  partitioned  into  a  region  in  which  all  messages  have  been  success¬ 
fully  transmitted,  and  a  region  that  has  been  returned  to  the  waiting 
interval. 

The  renewal  state  Incorporated  into  Gallager's  algorithm  is 
significant.  This  implies  that  the  system  repeatedly  comes  to  a  point 
in  the  algorithm  where  the  channel  history  is  independent  of  the 
statistics  of  any  interval  that  will  be  transmitted  in  the  future. 

This  is  a  feature  of  the  algorithm  that  we  will  try  to  maintain  when 
errors  are  introduced  into  the  system.  The  Gallager  procedure  is  not 
based  on  storing  the  result  of  slot  transmissions  from  the  beginning  of 
time,  but  only,  at  worst,  from  the  most  recent  renewal  state.  This 

makes  implementation  of  the  Gallager  algorithm  very  manageable. 

1  266 

Gallager  has  shown  that  the  optimal  value  of  £.q  is  — ^ —  and 
that  a  capacity  of  .48711  can  be  achieved  with  his  algorithm.  We 
note  that  this  is  a  FCFS  algorithm  in  that  no  message  is  successfully 
transmitted  before  a  message  with  an  earlier  generation  time. 

The  Increase  in  capacity  of  the  Gallager  algorithm  over  the 
Capetanakis  tree  algorithm  can  be  attributed  to  two  factors.  First 
of  all,  when  no  messages  are  found  in  the  first  half  of  a  conflict 
interval  Gallager  defines  the  second  half  of  the  interval  as  a  new 
conflict  interval  and  then  transmits  only  its  first  half  on  the  next 
slot.  In  the  analogous  situation,  Capetanakis  transmits  the  entire 
second  half  of  an  original  conflict  branch.  This  is  highly  inefficient 
(it  will  always  create  a  conflict). 


Secondly,  Capetanakis'  root  node  is  constrained  to  have  an 
integer  number  of  branches.  In  the  Gallager  algorithm  the  maximum 
transmission  interval  length,  can  be  optimized  over  a  continuous  range 
of  values.  This  is  also  a  source  of  superior  performance. 

2.2  The  Markov  Process  Model 

As  is  perhaps  implied  by  my  description  of  the  Gallager  algorithm, 
this  procedure  can  be  represented  as  a  Markov  process  with  a  countably 
infinite  number  of  states.  After  any  slot  we  are  faced  with  processing 
one  of  three  types  of  intervals:  an  interval  of  length  H  of  which  we 
have  no  knowledge  (after  a  renewal  state);  an  interval  which  we  know 
to  contain  more  than  one  message  (after  a  conflict);  or  an  interval 
which  we  know  to  contain  one  or  more  messages  (after  a  single  message 
has  been  found  in  the  first  half  of  a  conflict  interval).  We  thus 
define  the  following  state  classes: 

S  (T)  :  The  renewal  state  in  which  an  interval 
o 

of  length  T  will  be  transmitted  during 
the  next  slot. 

S^(x)  :  An  interval  of  length  x  contains  >_  1 

messages.  The  entire  interval  is  transmitted 
during  the  next  slot. 

S^(x)  :  An  interval  of  length  x  contains  2 

messages.  The  first  half  of  the  interval  will 
be  transmitted  during  the  next  slot. 

The  time  variables  are  normalized  here  by  multiplying  by  A.  Thus 
T  *  \iQ  and  is  in  units  of  packets.  Similarly,  x  is  in  units  of 
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packets  and  must  be  of  the  form  (T ) (2  )  where  i  is  a  nonpositive 

integer  (hence  the  states  are  countably  infinite).  1  will  continue 

to  use  this  notation  to  avoid  explicitly  carrying  X.  The  "length"  x 

is,  in  reality,  an  interval  where  the  expected  number  of  messages  is  x. 

Figure  2  shows  the  state  transition  diagram  for  the  Harkov 

process  model  of  the  Gallager  algorithm.  We  distinguish  between  the 

two  S^Oj)  states  stemming  from  a  S£(x)  state.  S0(-)j  indicates  that 

the  interval  being  considered  is  the  first  half  of  the  last  conflict 

interval  (the  second  half  has  been  returned  to  the  waiting  interval), 
x 

*>2^2  indicates  that  the  interval  being  considered  is  the  second  half 
of  the  last  conflict  interval  (the  first  half  has  been  successfully 
processed).  The  policy  for  the  two  states  is  identical. 

2.3  A  More  General  Algorithm 


Mosely  1 3 J  has  examined  the  Gallager  algorithm  with  an  eye  towards 
generalizing  it  in  two  respects.  The  first  issue  is  that  Gallager 
arbitrarily  splits  conflict  intervals  in  half,  when  the  optimum  dividing 
point  would  be  a  function  of  the  interval  length  and  X.  Secondly,  in 
a  S^(x)  state  Gallager  sends  the  entire  interval.  In  fact,  it  may  be 
more  efficient  to  send  only  a  portion  of  the  interval,  or  an  interval 
of  length  greater  than  x. 

Mosely's  analysis  determined  that  Gallager's  S^(x)  policy  was 
optimum.  She  also  found  that  splitting  conflict  intervals  into  equal 
parts  was  a  very  good  approximation  to  the  optimum  policy.  Whatever 
loss  of  optimality  that  is  Incurred  by  continuing  to  work  from  Gallager's 
algorithm  is  more  than  offset  by  the  ease  of  analysis  and  implementation 
it  provides. 


Markov  State  Transition  Diagram  for  Gallager's  Algorithm 


III.  THE  NOISY  CHANNEL  ALGORITHM 


3.1  Introducing  System  Errors 

Thus  far  we  have  assumed  that  all  sources  have  received  perfect 
information  while  monitoring  channel  activity.  In  other  words,  if  a 
source  detects  a  conflict  during  a  slot  transmission,  more  than  one 
message  was  sent  (with  probability  1).  In  this  chapter  a  limited 
class  of  probabilistic  detection  errors  is  introduced  into  the  system. 

We  then  examine  how  these  errors  influence  the  functioning  of  the 
Gallager  algorithm.  The  case  now  exists  where  updated  state  policy 
will  be  based  on  information  that  inaccurately  represents  the  statistics 
of  the  intervals  being  processed.  This  creates  a  number  of  system 
error  states  that  must  be  characterized. 

After  the  system  error  states  have  been  incorporated  into  the 
Markov  process  model  the  Gallager  algorithm's  ability  to  process  these 
states  is  examined.  It  becomes  apparent  that  the  algorithm  no  longer 
represents  a  stable  channel  allocation  scheme.  A  number  of  modifications 
are  then  made  to  the  Gallager  procedure  to  formulate  the  noisy  channel 
conflict  resolution  algorithm.  This  new  algorithm  retains  some 
important  features  of  the  Gallager  algorithm  and  possesses  a  number 
of  error  correcting  mechanisms. 

In  modeling  detection  errors  we  begin  by  assuming  that  when  an 
error  is  made,  all  sources  and  the  receiver  make  the  same  error.  No 
isolated  source  errors  occur,  only  system-wide  errors.  This  is  realistic 
if  the  sources  and  receiver  are  relatively  close  together,  or  if  the 
sources  receive  feedback  information  from  the  central  facility.  This 
assumption  allows  us  to  continue  to  model  the  system  as  a  Markov  process 


where  there  is  a  single,  system-wide  state  at  any  slot.  Sources  making 
errors  independently  would  create  a  situation  where  different  message 
generation  intervals  would  be  transmitted  during  the  same  slot.  This 
would  seriously  threaten  the  ability  of  a  transmission  interval 
algorithm  to  properly  resolve  conflicts.  The  analytical  approach  in 
dealing  with  isolated  errors  would  be  to  look  for  means  by  which  a 
source  that  has  made  an  error  could,  by  continued  monitoring  of  the 
channel,  return  to  the  correct  system-wide  state. 

There  are  six  types  of  detection  errors.  They  are 

PqI  =  Pr (  1  detected  |0  sent) 

Pq2  =  ?r  (  >_  2  detected  |0  sent) 

P^0  *  Pr  (  0  detected  |l  sent) 

P^2  *  Pr  (  2  detected  |l  sent) 

P2Q  =  Pr  (  0  detected  |  2  sent) 

P22  =  Pr  (  1  detected  |  2  sent) 

We  will  use  this  P. .  notation  to  indicate  both  a  type  of  error  and  a 
rj 

probability. 

The  Noisy  Channel  Algorithm  is  effective  if  any  of  the  above 
errors  occur.  However  its  analysis  appears  extremely  complicated 
unless  we  assume  PQ1  =  P^^  «  P2^  *  P2Q  =  0.  This  is  now  justified. 

With  all  six  error  probabilities  taking  on  non-zero  values  a 
great  number  of  error  states  are  introduced  into  the  Markov  representation. 
The  analysis  of  any  conflict  resolution  algorithm  rapidly  grows  in 
complexity  as  a  result.  It  can  also  be  argued  that  Pq2  and  P^  are  two 
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of  the  more  likely  errors  to  occur.  It  is  unlikely  that  noise  would  be 
interpreted  as  a  packet  with  the  necessary  protocol  and  parity  bits 
(P0i>.  It  is  also  unlikely,  from  a  signal  to  noise  ratio  argument, 
that  a  message  would  be  undetected  (P2q  or  ^10^ • 

The  underlying  sensitivity  to  channel  noise  can  be  determined  by 
measuring  the  system  capacity  with  varying  PQ2  and  P^.  Tnis  is  the 
important  issue. 

A  number  of  nice  properties  follow  from  the  P^.,,  P  only  error 
model.  These  include: 

1.  An  error  can  only  occur  when  a  conflict  is  detected  or, 
conversely,  a  detection  of  0  or  1  message  is  always 
correct. 

2.  All  errors  eventually  lead  to  "not  finding  enough  messages". 
These  can  be  discovered  as  processing  continues. 

3.  There  is  no  error  cancellation.  There  is,  for  instance, 
no  way  a  P^q  error  can  negate  a  Pq2  error. 

These  properties  contribute  toward  making  system  error  detection  and 
correction  a  workable  problem. 

There  are  now  three  error  state  classes  in  addition  to  the 
original  three  Gallager  state  classes: 

S2o (x)  :  An  interval  of  length  x  is  thought  to  contain  more 

than  one  message  but  contains  none.  This  is  created 

from  the  previous  slot  by  a  P  error  occuring  during 

transmission  of  the  x  interval,  or  by  no  message 

being  detected  (correctly)  in  the  first  half  of  a 
S2q(2x)  interval. 
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S21(x)  :  An  interval  of  length  x  is  thought  to  contain  more 
than  one  message  but  contains  exactly  one  message. 
This  arises  from  a  P^2  error  occuring  during  trans¬ 
mission  of  the  x  interval,  or  by  no  message  being 
found  in  the  first  half  of  a  S^(2x)  interval. 

S^q(x)  :  An  interval  of  length  x  is  thought  to  contain  one 

or  more  messages  but  contains  none.  This  arises  from 
a  single  message  being  found  in  the  first  half  of  a 
S2^(2x)  interval. 

It  is  important  to  realize  that  for  any  interval  conflict  resolution 
scheme  the  policy  at  each  state  is  only  dependent  on  the  number  of 
messages  thought  to  be  contained  in  an  interval.  In  other  words, 
the  policy  will  be  the  same  for  the  states  S21(x),  S2^(x),  and  S2(x). 
The  same  is  true  for  S^fx)  and  S^(x)  . 

Gallager's  algorithm  is  not  stable  for  the  ?q2’^12  c^ass 
errors.  If  a  Pq2  error  is  made  the  S2^(x)  state  is  entered.  The 
system  then  proceeds  to  keep  splitting  the  interval  in  search  of 
messages  that  don't  exist.  Unless  another  error  occurs,  a  new 
conflict  interval,  equivalent  to  the  most  recent  half  of  the  previous 
conflict  interval,  is  defined  after  each  slot.  If  another  P^2  error 
does  occur  during  processing,  the  system  still  progresses  to  a 
S2q(x)  state.  In  either  case,  no  packets  will  again  be  successfully 
transmitted. 

It  can  be  argued  that  if  the  Pq^  error  had  been  included  in  the 
model  this  type  of  instability  would  not  necessarily  occur.  The 
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important  point,  however,  is  that  the  Gallager  algorithm  has  no  way 
of  realizing  it  is  in  a  S^qCx)  state.  This  is  inefficient  and  some 
alternate  policy  at  S2  must  be  developed.  A  policy  must  also  be 
defined  for  the  case  where  a  transmission  at  results  in  no  message 
being  found. 

A  more  subtle  problem  also  arises  in  the  context  of  using  the 
Gallager  algorithm  to  process  detection  errors.  Suppose  the  system 
is  in  a  S^Cx)  state  and  a  error  occurs.  The  algorithm  states  that 
when  a  conflict  is  detected  in  the  first  half  of  a  conflict  interval, 
the  second  half  is  returned  to  the  waiting  interval.  This  is  based  on 
the  assumption  that  the  interval  being  returned  has  a  Poisson  message 
distribution.  In  the  case  indicated,  however,  the  interval  being 
returned  contains  two  or  more  messages.  The  claim  that  every  portion 
of  the  waiting  interval  has  an  identical  message  distribution  is  no 
longer  valid.  The  existence  of  a  renewal  state,  in  the  sense  that 
was  discussed  earlier,  is  also  lost.  We  no  longer  have  regular 
instances  where  the  channel  history  is  divorced  from  future  processing. 

The  Noisy  Channel  Algorithm  closely  follows  the  strategy  of  the 
Gallager  algorithm.  The  algorithm  has  been (modified ,  however,  to  include 
two  new  procedures.  One  of  these  insures  that  any  interval  being  re¬ 
turned  to  waiting  has  a  Poisson  distribution  of  messages.  This 
preserves  the  renewal  state.  The  other  new  policy  provides  a  means 
for  detecting  from  an  S2q(x)  state  that  an  error  has  been  made.  This 
restores  system  stability. 

3.2  The  Noisy  Channel  Algorithm 

In  this  algorithm  a  stack  mechanism  is  proposed  to  avoid  returning 
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"non-Poisson  intervals”  to  the  waiting  interval.  A  stack  entry  will 
consist  of  three  pieces  of  information:  the  length  of  an  interval,  a  lag, 
and  the  number  of  messages  that  have  been  successfully  transmitted  since 
that  stack  entry  was  created.  The  lag  must  be  incremented  after  each 
slot  and  the  message  counter  must  be  incremented  with  each  successful 
transmission.  All  intervals  that  would  normally  have  been  released  to 
the  waiting  interval  are  placed  on  the  stack.  Thus  an  entry  will  be 
added  to  the  top  of  the  stack  when  we  are  in  a  $2(x),  S^^Cx),  or  S £g(x) 
state  and  a  conflict  is  detected  upon  transmission.  All  stack  entries 
are  disjoint  and  the  top  stack  entry  represents  the  stack  interval 
furthest  removed  from  real  time. 

When  the  message  counter  for  a  stack  entry  reaches  two  that 
interval  is  released  to  the  waiting  interval  and  the  stack  entry  is 
eliminated.  The  two  successful  transmissions  imply  that  the  system 
was  in  fact  in  S^,  and  an  error  induced  conflict  did  not  occur  when 
that  stack  entry  was  created.  We  verify  that  an  interval  has  a  Poisson 
message  distribution  before  it  is  released  to  the  waiting  interval. 

Two  successful  transmissions,  as  opposed  to  a  detected  conflict,  are 
necessary  for  this  verification.  Stack  entries  will  be  eliminated  from 
bottom  to  top,  or  equivalently,  from  right  to  left  on  the  time  axis. 

The  stack  is  used  to  determine  the  next  state  when  an  error  is 
detected.  If  we  come  to  a  point  in  processing  where  we  are  "a  message 
short"  (a  transmission  yields  no  message)  then  we  take  the  top  entry 
off  the  stack  and  process  it  as  if  it  were  Poisson  conditioned  on 
having  1  messages .This  is  as  opposed  to  normally  processing  the 
interval  as  if  it  were  Poisson  conditioned  on  having  >  0  messages. 
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If  Che  stack  is  empty  it  indicates  that  the  most  recent  detected  conflict 
was  a  error,  and  the  system  returns  to  S^.  S2J  errors  are  corrected 
in  the  same  manner. 

Similarly,  if  we  discover  the  system  is  in  a  $2^  state,  we  take 
the  top  entry  off  the  stack  and  follow  S2  policy.  If  the  stack  is  empty 
the  last  detected  conflict  was  a  error  and  we  return  to  Sq. 

If  the  system  is  in  a  regular  state  (S^.S^,  or  S2)  it  indicates 
that  all  processing  errors  have  thus  far  been  corrected.  In  the 
renewal  state,  S^,  the  stack  will  be  empty.  The  or  S2  state,  however, 
does  not  indicate  an  empty  stack. 

The  next  question  that  arises  is  how  does  the  algorithm  "discover" 
that  it  is  in  a  S2Q  or  $2^  state?  The  logical  approach  to  this  problem 
is  that  at  some  point  we  want  to  stop  splitting  a  conflict  interval 
and  transmit  the  entire  interval.  If  the  conflict  is  confirmed  then  the 
splitting  procedure  is  continued.  If  an  idle  channel  or  a  successful 
transmission  is  detected  then  an  error  must  be  corrected  through  use  of 


the  stack.  This  is  the  strategy  of  the  Noisy  Channel  Algorithm. 


We  define 


_  -NCTOFF 
x  ■  T  •  2 


where  NCTOFF  is  a  nonnegative  integer.  We  then  specify  a  new  policy  for 


states  S2q(x),  S21(x),  and  S2(x)  when  x  <_  xc. 


For  x  x^  the  system  instantaneously  moves  from  these  states 

(before  the  next  slot  transmission)  to  a  special  processing  state 
*  * 

S^(x).  The  policy  at  S^(x)  calls  for  the  transmission  of  the  entire 


interval.  If  0  or  1  message  is  detected  then  a  new  state  is  determined 
using  the  stack.  An  error  has  been  detected.  If  a  conflict  is  detected 


during  S^(x)  transmission  then  the  system  moves  to  another  special 

£  A 

processing  state  -  S^Cx).  The  policies  at  S2(x)  and  S2(x)  are  the  same, 
except  that  the  S2(x)  policy  disregards  the  x£  cutoff.  Therefore  the 
conflict  interval  will  be  split  once.  If  a  conflict  is  detected  the 

X  *  X 

system  will  go  to  state  which  will  in  turn  go  to  ^(-jr).  ^us 

the  system  going  to  $2(x)  allows  one  interval  division  where  the 
policy  at  S2<x)  would  not. 

If  there  are  in  fact  two  messages  contained  in  an  interval  of 
length  x  then  this  policy  change  is  not  efficient.  The  first 

half  of  every  new  conflict  interval  will  be  transmitted  for  two 
consecutive  slots  until  the  conflict  is  resolved.  A  sort  of  oscillation 
between  and  S2  policy  will  take  place.  An  optimum  xc  must  therefore 
be  determined  so  that,  based  on  an<1  P22*  t*ie  risk  transmitting 
intervals  thought  to  contain  ^  2  messages  is  justified. 

This  policy  for  x  xc  can  be  viewed  as  a  member  of  a  general 
class  of  algorithms.  We  could  follow  a  procedure  of  sending  an  entire 
S„(x)  interval  for  x  <  x  and  then  allow  for  n  divisions  before  a 
conflict  interval  would  again  be  transmitted  in  its  entirety.  In  the 
algorithm  above  n  =  1.  A  larger  n  would  make  sense  for  a  system  with 
small  detection  error  probabilities,  but  it  is  doubtful  the  capacity  would 
be  significantly  higher  than  the  for  the  n  *  1  policy,  with  NCTOFF 
and  T  optimized. 

As  an  example  of  the  Noisy  Channel  Algorithm,  in  Figure  3  we 
have  the  same  message  processing  problem  that  was  used  in  Figure  1  to 
illustrate  the  Gallager  Algorithm.  A  error  has  occurred  in 
figure  3.2b,  however,  at  a  state.  We  assume  that  NCTOFF  is  set 


at  three  for  the  algorithm.  At  figure  3.2c  the  interval  A  has  been 

U 

put  on  the  stack  and  the  system  is  following  the  policy  of  the  ®20 ( — 4^) 
state.  Figures  3.2c  and  3. 2d  show  succesive  conflict  interval  divisions 

_3 

yielding  no  messages.  In  figure  3.2e  the  cutoff  point  (2  )  (X8.  ),  has 

l  ° 

been  reached  and  the  entire  interval  of  length  — ^  slots  is  transmitted. 

When  this  creates  an  idle  channel  the  interval  A  is  taken  off  the 

i  \ 

stack  and  the  system  is  in  the  state  ^ — )  in  figure  3.2f.  Processing 

then  continues  as  in  the  Gallager  Algorithm.  We  have  introduced  only 
one  error  into  this  example.  Another  error  could  have  occurred  at  any 
point  in  the  algorithmic  processing  where  0  or  1  messages  was  transmitted. 
This  would  have  created  the  necessity  for  a  greater  number  of  slots  to 
successfully  transmit  the  messages.  An  error  is  detected  and  then  we 
"backtrack"  by  taking  entries  off  the  stack,  to  the  point  where  the  error 
was  made. 

The  noisy  channel  conflict  resolution  algorithm  is  explicitly 

defined  in  Appendix  A  by  the  Markov  state  transition  policies.  Note 

that  while  Gallager 's  state  is  depicted  in  this  model,  it  now  has 

a  different  associated  policy  that  involves  interaction  with  the  stack 
* 

and  .  A  policy  is  also  defined  for  the  case  where  a  transmission 
from  a  state  results  in  no  messages  being  detected.  The  policy  at 
Sq  is  unchanged. 

The  state  of  the  system  at  any  slot  is  now  also  a  function  of  the 
stack  entries  at  that  point.  The  system  state  has  an  additional 

dimension.  A  reference  to  a  system  state,  like  is  really  a 

reference  to  a  class  of  states.  The  states  in  a  class  are  differentiated 
by  their  interval  length  and  the  stack  configuration. 


The  Noisy  Channel  Algorithm  corrects  all  and  errors 
and  is  stable  for  this  class  of  errors.  Certain  features  of  the 
Gallager  Algorithm  have  been  preserved.  Namely,  the  Noisy  Channel 
Algorithm  is  FCFS  and  has  a  renewal  state. 

This  algorithm  must  now  be  optimized  over  two  parameters:  T  and 
Xc(or  equivalently,  NCTOFF) .  A  quantitative  discussion  of  the  algorithm 
and  this  optimization  is  presented  in  the  next  chapter. 

IV.  THE  CAPACITY  CALCULATION 
4.1  The  General  Approach 

An  approach  to  calculating  the  capacity  of  the  noisy  channel 
system  was  hinted  at  in  chapter  2.  Just  as  in  the  Gallager  algorithm, 
the  expected  change  in  the  system  lag  between  successive  entries  in 
state  SQ  must  not  be  positive  in  the  Noisy  Channel  Algorithm.  A 
stable  system  cannot,  on  the  average,  be  successfully  transmitting 
messages  at  a  rate  slower  than  the  message  generation  rate.  Equivalently, 
the  waiting  interval  cannot  have  a  positive  expected  growth  rate. 

In  this  chapter  we  derive  an  expression  for' the  expected  change  in  the 
system  lag  between  successive  renewal  states.  This  is  a  linear,  monoto- 

ft 

nically  decreasing  function  in  1/X,  and  therefore  the  X  for  which  the 
expected  lag  change  equals  zero  is  the  system  capacity. 

Gallager* s  algorithm  involved  defining  a  single  policy  in  terms 
of  T*-  the  normalized  transmission  interval  that  yielded  the  highest 
stable  throughput  (the  capacity).  With  the  Noisy  Channel  Algorithm  the 
optimum  policy  is  a  function  of  P^  and  P^* 


In  addition,  there  are 
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two  parameters,  T  and  NCTOFF,  that  must  be  jointly  optimized  for  each 
set  of  error  values.  We  will  want  to  examine  the  manner  in  which  the 


optimum  policy,  as  well  as  the  capacity,  varies  with  different  system 
error  probabilities. 

4.2  The  Capacity  Calculation 

We  begin  by  defining  the  following  random  variable: 

L(Sq)  :  The  change  in  the  system  lag  from  a  renewal  state  slot 
to  the  next  renewal  state  slot. 

The  expected  lag  change,  E(L(SQ)),  has  three  components: 

E(L(Sq))  *  Etnumber  of  slots  until  the  next  renewal  state] 

+  E[number  of  slots  returned  to  the  waiting  interval] 

'  I  <*•» 

is  the  initial  Sq  transmission  interval  length  -  in  units  of  slots. 
This  is  the  time  that  has  been  taken  from  the  waiting  interval. 

We  next  decompose  L(Sq)  into  five  terms.  These  correspond  to 
L(Sq)  conditioned  on  the  five  possible  sent-detected  message  combinations 
that  could  occur  during  the  first  transmission  interval.  We  define  these 
random  variables  as  follows: 

Lqq  :  The  system  lag  change  between  successive  renewal  states 
given  that  the  first  transmission  results  in  0  messages 
sent  and  0  messages  detected. 

Lq2  :  The  system  lag  change  between  successive  renewal  states 
given  that  the  first  transmission  results  in  0  messages 
sent  and  >  2  messages  detected. 
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states  given  that  the  first  transmission  results  in 
1  message  sent  and  1  message  detected. 

1*12  !  The  system  lag  change  bet  »een  successive  renewal  states 
given  that  the  first  transmission  results  in  1  message 
sent  and  2  messages  detected. 

^22  :  system  J-a8  change  between  successive  renewal  states 

given  that  the  first  transmission  results  in  >  2  messages 

and  2  messages  detected. 

It  follows  that, 

E(L(SQ))  -  Pr(0  in  T)  [PQ2  E(L02)  +  (1  -  PQ2)  E(LQq) ] 

+  Pr  (1  in  T)  [P12  E(L12)  +  (1  -  P^)  E(Ln)] 

+  Pr(>  2  in  T)  E(L22)  (4.2) 

Clearly  E(L^j)  can  be  thought  of  as  having  the  same  three  types 
of  contributing  terms  as  E(L(Sq)).  The  only  instance  where  the  possibility 
of  returning  time  to  the  waiting  interval  exists  is  when  >2  messages  are 

contained  in  the  initial  interval.  Thus  only  E[L22J  will  have  a  non¬ 

zero  return  term. 

E(Lqo)  and  E(L22)  are  straight-forward  to  calculate.  These 
correspond  to  the  case  where  all  the  messages  in  an  interval  of 
normalized  length  T  are  successfully  transmitted  in  one  slot.  There¬ 
fore 


E'V 

■ 

(4.3a) 

E'Llll 

-  l-f 

(4.3b) 
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In  considering  the  other  three  terms  we  define  the  following 
functions : 

t„^(x)  :  The  expected  number  of  slots,  assuming  an  initially 

empty  stack,  before  a  renewal  state  is  reached  from 

the  state  S2q(x). 

t^^Cx)  :  The  expected  number  of  slots,  assuming  an  initially 

empty  stack,  before  a  renewal  state  is  reached  from 

the  state  S^^Cx). 

t22(x)  •'  The  expected  number  of  slots,  assuming  an  initially 

empty  stack,  before  a  renewal  state  is  reached  from 

the  state  S22(x). 

R(x)  :  The  expected  interval  length  returned  to  waiting  from 

an  interval  of  length  x  known  to  contain  _>  2  messages. 


In  t20(x),  t21(x)  and  t22(x)  the  interval  length  x  has  again 
been  normalized.  It  represents  a  number  of  slots  multiplied  by  the  arrival 
rate.  These  intermediate  functions  are  therefore  independent  of  A.  The 
empty  stack  assumption  in  these  definitions  implies  that  exactly  one 
conflict  has  been  detected  since  the  previous  renewal  state. 

From  these  definitions  it  follows: 


E(L02>  ’ 

(1  + 

t20(T» 

T 

'  X 

(4.4a) 

e<l12)  - 

(1  + 

t21(T)) 

T 

"  X 

(4.4b) 

e(l22>  ■ 

(1  + 

t22(T)) 

_  T  -  R(T) 

X 

(4.4c) 

In  each  equation  here  the  1  corresponds  to  the  initial  slot  that 
sends  the  system  into  one  of  the  three  S2~type  states. 

Once  t2Q»  t2i»  and  R  have  been  determined  the  capacity 
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calculation  is  essentially  done. 

In  computing  t2g(x),  tj^x)  and  t22(x)  there  are  three  important 
factors.  First  of  all,  these  functions  will  be  defined  recursively  in 
terms  of  the  expected  times  to  process  intervals  of  length  ^  (created 
by  the  splitting  of  a  conflict  interval).  In  each  case  we  look  at  the 
policy,  the  different  possible  message  distributions,  and  the  different 
possible  sent-detected  pairs  (based  on  PQ2  and  P^)-  Each  oi  the  functions 
is  in  the  form  of  a  constant  term  plus  probabilistic  weighting  factors 
multiplied  by  t 2q «  *21^’  an<*  t22^'  Secondly,  processing  policy 
changes  for  x  <_  and  we  would  expect  the  recursion  formulae  to  also 
change  at  that  point.  This  is  where  NCTOFF  enters  the  capacity  calculation. 
Thirdly,  a  more  subtle  point,  although  we  have  defined  these  functions 
assuming  an  empty  stack,  all  stack  manipulations  are  implicitly  accounted 
for  in  this  recursive  construction.  This  will  be  explained  in  more  detail 
later. 


We  first  consider  t~n(x)  for  the  case  where  x  <  x  .  Both  the 
zu  —  c 

message  distribution  (there  are  no  messages)  and  the  policy  followed 

are  independent  of  x.  t2Q(x)  for  x  £  x^  is  therefore  dependent  on  only 

* 

Pq2  and  P^2‘  Call  it  t2Q  •  We  now  consider  exactly  what  happens  in  the 
Noisy  Channel  Algorithm  when  we  have  an  interval,  of  length  x  _<  xc  where 
we  think  there  are  2  messages,  but  there  are  none.  The  entire  interval 
is  first  transmitted  and  this  will  bring  the  system  to  unless  a  Pq2 

error  is  made.  If  a  ?  ^  error  is  made  the  first  half  of  the  interval 
is  sent  next.  This  will  completely  process  (— )j  unless  another  P^j 
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error  forces  the  system  into  the  state.  In  any  case,  the 

system  will  eventually  determine  that  there  are  no  messages  in  (^)^, 
and  the  system  will  enter  S2q(^)2*  It  follows  from  this  discussion: 


1  +  P02  11  +  P02  t20  +  t201 
1  + 


1  -  P  -  P 

02  02 


(4.5a) 


(4.5b) 


Thus  the  algorithm  remains  stable  only  for  P^  < 
* 

the  point  at  which  t,_  goes  to  infinity. 


-1  + 


For  t^Q(x)  where  x  >  x£  the  policy  followed  is  no  longer  independent 
of  x.  The  closer  x  is  to  x^,  the  closer  we  are  to  a  policy  change. 

This  is  important  because  the  only  way  we  discover  the  system  is  in  a 
S2o(*)  state  is  to  transmit  the  entire  interval. 

In  processing  a  S2Q(x)  interval  the  first  half  of  the  interval  is 
transmitted.  If  this  results  in  a  Pq2  error  then  we  have  two  S2q(^) 
intervals  to  process.  (The  second  one  is  removed  from  the  stack  and 
processed  when  the  first  one  is  recognized  as  being  an  error.  This 
illustrates  the  manner  in  which  the  handling  of  the  stack  is  implicitly 
accounted  for  in  this  development).  If  no  error  occurs,  there  is  only 
S20^2  to  Proceas-  Therefore, 

t20(x)  “  1  +  (1  -  PQ2)  t2Q(-j)  +  2Pq2  t^Oj),  X  >  xc  (4.6a) 


t20(x)  *  1  +  (1  +  PQ2)  t2Q(-j)  , 


x  >  *c  (4.6b) 


We  now  have  a  procedure  for  calculating  t20^’  We  8tart  with 
t20^xc^  "  c20  an<*  use  t*'e  recurs^on  formula  of  equation  4.6b  NCTOFF 
times  to  reach  t2(j(T). 


‘21 W 

As  with  t2g(x),  t2^(x)  is  independent  of  x  for  x  £xc»  If 
there  is  one  message  in  an  interval  of  length  x,  it  is  in  either  half 
of  the  interval  with  probability  1/2. 

"k 

In  the  derivation  of  t-.  and  t., (x)  x  >  x  ,  we  use  the  fact 

*  1  Cl  C 

that  the  expected  time  to  reach  Sg(x)  from  an  empty  stack  S1Q(x) 
state  is  1  +  PQ2  t20(x). 

In  ^^(x)  w^ere  x  <  xc  we  first  send  the  entire  interval  which 
will  bring  the  system  to  Sg  unless  a  error  is  made.  If  an  error  is 
made,  the  first  half  of  the  interval,  is  transmitted  next.  We  then 

consider  two  equiprobable  cases:  contains  the  message  or  it  doesn't. 

If  the  message  is  in  (y)  ^  then  we  will  reach  a  S21^1  state  with  proba¬ 
bility  P^2»  and  we  will  eventually  reach  a  s^q(^)2  state  after  (^)^  has 
been  processed.  If  the  message  is  in  (^>2  we  will  reach  S2q(^)^  w*th 
probability  PQ2,  and  once  is  processed  we  will  reach  s21(^)2- 

Thus, 


-  1  +  p!2  [i  =  2  (Pi2  C21  +  1  +  P02  t20) 

+  2  (Pg2  t2Q  +  t21)] 

.  3  * 

1  +  —  p  +p  p  t 

2  12  *12  "p2  C20 

x=  - - - —  .  ■  >  - 

2 

P  +  P 
,  r12  12 


(4.7a) 


(4.7b) 


Thus  the  algorithm  is  stable  for  all  P^  <  1. 


tji^)  for  x  >  xc  is  no  longer  independent  of  x.  We  note  that 
t2j(x)  for  x  >  xc  can  be  related  to  t2^>  The  policies  are  identical  except 
that  there  is  no  initial  transmission  of  the  entire  interval.  C2i^x^  now 
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begins  by  transmitting  Also  t22(^)  must  be  substituted  for  t21> 

and  t2Q(|)  must  be  substituted  for  t*Q  in  the  expression.  This  same 
type  of  transformation  also  relates  t2Q  and  t2Q(x)  for  x  >  xc>  as  well 
as  t22(x)  x  _<  xc  and  t22(x)  x  >  xc- 

t2i^x)  =  \  l1  +  Po2  t20^  +  t21(r2.^ 

+  I  [1  +  P12  t21(!)  +  1  +  P02  t20(!)1  *  X  >  Xc  (4*8a) 


t21^X^  =  2  +  Pn,)  + 


1  +  P 


i-2  t  (2-) 

2  c2lV  ' 


x  >  x  (4.8b) 
c 


02  20  2 

Thus  we  have  a  straightforward  method  to  recursively  compute 


t22(x) 

We  now  consider  where  x  ^_xc*  an^  realize  that  this 

function  is  not  independent  of  x.  The  distribution  of  >_  2  messages 

in  an  interval  is  dependent  on  its  length.  For  instance,  as  x  approaches 

zero  we  would  expect  that  there  are  exactly  two  messages  in  the  interval. 

For  large  x  there  would  be  a  much  greater  probability  of  having  more 

than  two  messages  in  the  interval.  t22(x),  f°r  X  —  Xc’  is  3  funct:LOn  of 

* 

t*  ,  and  t00($).  Through  reasoning  similar  to  that  used  with  t2Q 
20  21  22  2 

and  t2*  it  follows  that: 

t22(x)  *  1  +  Pr(0  in  (^)]j:L  2  in  x)  I1  +  P02  fc20  +  C22^2^ 

+  Pr  (1  in  (^)  xi  >_  2  in  x)  • 

[2  +  Pr(>2  in  (|>2I  >1  in(|)2>  •  t22<|)  +  (1  +  Pr  (1  in  (f )  2I 

in  (|)2))  Pl2  t21]  +  Pr(>2  in  (f^ll2  i«  X)U  +  t22(f)]* 

x  <  x  (4.9) 


h  Unr  i  iminrt*tnrr 
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We  see  that  after  the  initial  transmission  (which  results  in  a 
conflict)  there  are  three  terms  that  arise  from  the  possible  number  of 
messages  in  the  first  half  of  the  interval.  If  (-j)^  contains  no  messages 

X  X 

then  S22^2  wi^  eventually  be  reached.  If  (j)^  contains  ^  2  messages 

X  X 

then  (^2  returned  to  waiting  and  only  S22  ^2*1  nee(*8  to  processed. 

X 

If  there  is  exactly  one  message  in  (-j)^  then  two  cases  are  considered: 

X  X 

(-j) 2  containing  exactly  one  message  and  containing  more  than  one 

message . 

t22(x)  for  x  >  xc  follows  directly  from  the  previous  derivation. 
There  is  no  initial  transmission  of  the  entire  interval,  and  t2o^  an<* 

X  k  k 

*"21^  are  substitutet*  f°r  fc20  an<*  C21  resPectively- 

t22(x)  »  Pr(0  in  (|)x|  >2  in  x)  [1  +  PQ2  t2Q(|)  +  t22(|)] 


+  Pr  (1  in  (|)1|  >2  in  x)  [2  +  Pr  (> 2  in  (|>2|  >1  in  (f>2)  * 

t22(|)  +  (1  +  Pr(l  in  (|)2|  >1  in  (|>2))- 

P12  C21(2)J 

+  Pr  (>  2  in  (|)x|  >2  in  x)  [1  +  t22(|)]  x  >  xc  (4.10) 


With  t2g(x)  and  t21(x)  we  had  an  initial  value  from  which  to 
k  k 

begin  recursion  (t2Q  or  t^).  This  is  no  longer  true  in  calculating 
t22(x).  The  approach  used  here  is  to  make  a  linear  approximation  to 
t22^J0  for  x£  x  -  0.  This  then  serves  as  an  initial  value.  One 
recursion  formula  is  used  to  calculate  t„_(x  ).  The  second  recursive 
procedure  is  then  used  to  get  to  t22(T).  The  linearization  procedure 
and  the  computational  accuracy  of  this  method  are  discussed  in  Appendix  B. 
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The  results  are  summarized  below: 


t22(x)  -  a  +  bx 

x  >  x  =  0 
c  — 

(4.11a) 

* 

P  t 

c  .  02  20 

3*5+  - 

2 

+  2P  t* 

12  21 

(4.11b) 

b  ■  18  119  +  p02 

C20  +  3P12  t21i 

(4.11c) 

We  now  turn  to  computing  R(T)  -  the  expected  interval  length 
returned  to  waiting.  This  is  a  function  of  the  message  distribution 
in  T,  and  is  independent  of  the  error  probabilities.  R(x)  will  be 
recursively  defined  in  terms  of  RCj).  If  there  are  less  than  two 
messages  in  (y)^  and  two  or  more  messages  in  (-^ )  2  t*ien  we  consider 

X  XX 

R ("2 )  2  *  ^  there  are  two  or  more  messages  in  (-j)  ^  then  *s  returned 

to  waiting  and  we  determine  R(^)^.  We  see  that: 

R(x)  -  Pr (0  in  (|)1|  >2  in  x)  R(|) 

+  Pr  (1  in  (|)J  >2  in  x)  Pr(>2  in  (|)  2 1  >1  in  (|)2>  R(|) 

+  Pr(>  2  in  (j)x\  >2  in  x)  (|  +  R(|))  (4.12) 

We  note  that  R(x)  is  independent  of  x .  The  same  recursion  formula 
is  used  for  any  value  of  x.  We  again  have  the  problem  of  obtaining  an 
initial  value  and  a  linearization  procedure  is  used  for  x  a  0.  The 
details  of  this  approximation  are  given  in  Appendix  B. 

R(x)  =  a  +  bx  ,  x  =  0  (4.13a) 

a  =  0  (4.13b) 
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b  =  1/6 


(4.13c) 


We  now  have  defined  explicit  means  for  calculating  » 

C2i(T) »  and  R^-  Xt  follows  that  EILq23»  e[L121.  and  eIL22^ 

can  be  found  from  equation  4.4.  Substituting  into  equation  4.2 
yields  a  linear  equation  in  1/X  for  E[L(Sq)].  Setting  E[L(Sq)]  =  0 
and  solving  for  X  yields  the  largest  stable  throughput  that  can  be 
achieved  for  the  values  of  T  and  NCTOFF  that  have  been  used.  Appendix  C 
contains  the  FORTRAN  programs  that  were  used  in  making  this  calculation, 
including  the  recursive  iterations  involved  in  determining  the  inter¬ 
mediate  values  that  have  been  discussed. 

V.  ANALYSIS  OF  NOISY  CHANNEL  ALGORITHM  PERFORMANCE 
5.1  Capacity  Degradation  Due  to  System  Errors 

In  this  chapter  the  performance  of  the  Noisy  Channel  Algorithm 

with  respect  to  and  P^  is  discussed.  In  addition,  a  few  observations 

*  * 

are  made  as  to  how  the  optimal  algorithm  parameters  (T  and  NCTOFF  ) 
behave  with  varying  system  error  probabilities. 

Figure  4  is  a  graph  of  the  algorithm  capacity  at  different  error 
levels.  P^  and  P^2  are  equal  and  denoted  by  P.  Note  that  the  parameters 
are  not  the  same  at  each  point  on  the  graph,  but  rather  each  point 
represents  the  achievable  throughput  when  the  optimal  stragegy  is  imple¬ 
mented  for  that  error  probability. 

It  follows  from  Figure  4  that  the  capacity  can  be  very  well 
approximated  by  a  line  over  the  region. 1  —  P  =  .5.  This  line  has  a  slope 


of  about  .58.  At  both  ends  of  the  error  range  the  capacity  falls  off 
at  a  faster  rate .  When  errors  are  first  introduced  into  the  system  (as 
we  move  away  from  P  =  0)  channel  performance  is  downgraded  more  rapidly 


than  when  the  error  probability  is  increased  the  same  amount  for  an 
already  noisy  channel  (say  P  =.3). 

The  system  performance  also  deteriorates  quickly  as  P  approaches  .62 

* 

This  is  the  value  at  which  t£Q  goes  to  infinity  and  the  system  there¬ 
fore  becomes  unstable  for  any  policy  and  any  non-zero  throughput.  This 
can  be  interpreted  as  the  point  where  there  are  too  many  error  induced 
"messages”.  The  system  cannot  recover  from  the  time  spent  correcting 
these  errors. 

* 

Thus  the  A  vs.  P  curve  exhibits  a  reasonable  and  understandable 
behavior  in  all  regions. 

It  should  be  noted  that  if  the  P^  =  P-^  constraint  is  removed 
P^2  can  take  on  any  value  less  than  1  (see  4.7b)  and  the  system  will 
still  be  able  to  achieve  a  non-zero  stable  throughput.  This  is  related 
to  the  fact  that  a  S2^^x)  interval  is,  in  general,  more  efficiently 
processed  than  a  S20  (x)  interval.  Once  a  message  is  successfully 
transmitted  in  the  first  half  of  a  S2^  interval  the  system  moves  to 
the  state  class.  This  then  transmits  the  entire  remaining  interval. 

The  algorithm  policy  at  an  S2Q  state,  however,  leads  to  repeated  divisions 

of  a  conflict  interval  that  contains  no  messages.  Even  below  the  cutoff 
length,  this  inefficiency  underlies  S2Q  vs.  S21  processing. 

P^2  going  to  1  can  simply  be  interpreted  as  the  system  never  being 
able  to  successfully  transmit  a  message.  Either  P^  going  to  1  or 
Pq2  going  .62  drives  the  system  capacity  to  zero. 

Table  1  indicates  how  the  policy  that  achieves  capacity  varies 

with  error  probability.  Again  P  =  P  =■  P^*  T  has  been  optimized  to 
the  nearest  .05.  A  few  important  characteristics  of  the  optimal 
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P“P02,P12 

NCTOFF* 

* 

T 

X* 

0.00 

CO 

1.25 

.487 

0.05 

5 

1.45 

.424 

0.10 

3 

1.45 

.383 

0.15 

2 

1.40 

.350 

0.20 

2 

1.55 

.323 

0.25 

1 

1.35 

.297 

0.30 

1 

1.45 

.273 

0.35 

1 

1.60 

.247 

0.40 

0 

1.30 

.219 

0.45 

0 

1.40 

.188 

0.50 

0 

1.60 

.151 

0.55 

0 

1.95 

.103 

0.60 

0 

2.90 

.035 

TABLE  1 


policy  can  be  extracted  from  this  table. 

* 

First  of  all,  NCTOFF  is  a  monotonicaliy  decreasing  function 

of  P.  As  the  probability  of  error  increases,  the  algorithm  takes 

fewer  conflict  interval  divisions  before  deciding  to  transmit  an  entire 

conflict  interval.  Secondly,  for  each  interval  of  P  having  constant 
*  * 

NCTOFF  ,  T  is  monotonicaliy  increasing  in  P.  Thus  for  a  given  optimum 
NCTOFF  we  take  a  larger  initial  transmission  interval  (in  normalized 
time)  as  the  chance  of  error  increases.  The  optimality  here  is  again 
related  to  the  fact  that  slots  lead  to  more  efficient  error  correction 
than  do  S2Q  slots.  This  becomes  more  important  as  P  increases,  and  we 
therefore  transmit  a  renewal  state  interval  with  a  greater  expected 
number  of  messages.  Offsetting  this  is  the  possibility  of,  if  T  is 
too  large,  transmitting  an  interval  that  has  too  high  a  probability 
of  containing  2  messages. 

* 

In  Table  1  we  see  that  when  NCTOFF  is  decremented  there  is  a 

*  * 

corresponding  discontinuous  decrease  in  T  .  T  then  increases  with  P 

★ 

again  until  the  next  drop  down  in  NCTOFF  . 

* 

Figure  5  shows  the  error  probabilities  at  which  jumps  in  NCTOFF 
occur.  As  P  approaches  0  the  jumps  occur  with  increasing  frequency. 

•k 

At  P  ■  0,  NCTOFF  must  go  to  infinity.  This  is  the  error  -  free  case 

and  we  would  never  want  to  transmit  a  conflict  interval  in  its  entirety. 

The  point  P  =  .40  is  also  significant  in  Figure  5.  This  is  the  error 

* 

probability  at  which  NCTOFF  goes  to  zero.  At  P  *  .40  the  probability 
of  error  is  great  enough  that  the  optimum  policy  dictates  a  second 
transmission  of  the  entire  interval  if  a  conflict  is  detected  during 
the  initial  transmission.  All  intervals  are  at  or  below  the  cutoff 
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point  and  are  handled  using  the  special  processing  states. 

Figure  6  is  a  graph  of  the  maximum  stable  throughput,  for 

PQ2  *  P^2  “  .5,  over  a  range  of  NCTOFF  and  T  values.  The  optimal  policy 

*  * 

for  these  error  probabilities  is  T  *1.6  and  NCTOFF  *  0.  We  see  that 
for  NCTOFF  *  0,  X  is  quite  flat  over  1.2  £  T  £  2.6,  and  the  choice  of 
the  value  of  T  is  not  critical  to  achieving  high  throughput.  We  also 
note  that  the  maximum  on  the  NCTOFF  *  1  curve  is  very  close  to  capacity. 
Taking  a  larger  transmission  interval  and  allowing  for  one  conflict 
interval  division  is  roughly  equivalent  in  performance  to  taking  a 
shorter  interval  and  immediately  following  cutoff  policy.  We  note  that 
as  NCTOFF  decreases,  the  maximum  throughput  point  on  each  NCTOFF  curve 
is  at  a  smaller  T  value.  NCTOFF  >  2  seems  to  significantly  downgrade 
the  performance  of  this  high  error  system.  We  see,  however,  that  the 
Noisy  Channel  Algorithm  gives  some  flexibility  in  choosing  parameters 
to  attain  near-capacity  throughput. 

5.2  Algorithm  Performance  for  Small  Error  Probabilities 

In  this  section  we  first  consider  the  degenerate  case  Pg2  =  Pj^  * 

This  corresponds  to  an  error-free  system  in  which  the  optimal 

algorithm  never  uses  the  special  processing  states. 

* 

As  was  explained  previously,  NCTOFF  goes  to  infinity  in  this  case. 

The  Noisy  Channel  Algorithm  therefore  reduces  to  the  Gallager  algorithm 

*  * 

in  this  instance.  T  *  1.266  and  A  =  .48711  -  concurring  with 
Gallager 's  results. 

In  practice,  capacity  can  be  very  nearly  achieved  with  NCTOFF  10 


Figure  7  shows  the  maximum  stable  throughput  for  different  algorithm 


Algorithm  Performance  at 


NCTOFF 


parameters .  \  is  approximately  constant  for  NCTOFF  “  10  and 

.9  <  T  <_  1.9.  We  also  observe  that  the  NCTOFF  “  5  and  NCTOFF  *  10 
curves  closely  follow  each  other.  So  again,  there  is  a  range  of 
effective  algorithm  parameters.  We  see  that  as  NCTOFF  approaches 
zero  the  inefficiency  of  transmitting  conflict  intervals  (which  will 
always  contain  >  2  messages)  begins  to  take  its  toll,  and  this  leads 
to  severely  downgraded  performance. 

Although  t^CH  an<^  l21^T^  are  n0t  rea3^  meaningful  in  the  case 
where  “  P^  *  can  be  shown  that  they  reduce  to  simple  closed- 

form  expressions.  From  the  derivations  in  Chapter  4,  for  P^  “ 


11 20  *  1 


t20(x)  -  1  +  t2o^2^ 


and  therefore,  1 2Q (T)  =  1  +  NCTOFF. 


Similarly, 


t2i  *  1 


t21(x>  -]+<§>  t21(|) 


and  therefore, 

NCTOFF- 1 

t21(T)  “  1  +  £  2_j 
j-0 

lim  t21(T)  ”  3’ 

NCTOFF-*® 

It  can  also  be  shown  from  limiting  state  Poisson  distributions 
(as  x^O)  and  the  special  processing  state  policies,  that  t22(0)  ■  5 
for  P  -  0.  This  is  as  predicted  by  the  linear  approximation  to  t22(x) 
given  in  appendix  B. 


$  . .  “ 
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Table  2  indicates  the  behavior  of  the  optimum  policy  and  capacity 

as  P  moves  away  from  0.  NCTOFF  =  10  was  the  greatest  value  considered, 

and  T  was  optimized  to  the  nearest  .005.  NCTOFF*  decreases  very  quickly 

as  P  is  increased  (from  00  to  3  as  P  moves  from  0  to  ,1).  X  is 

approximately  linear  over  0  _<  P  £.1.  A  good  rule  of  thumb  in  this 
* 

region  is:  X  =  .487  -  P  . 


5. 3  Some  Interesting  Algorithm  Characteristics 

A  few  other  interesting  points  can  be  made  from  studying  data 
related  to  the  functioning  of  the  Noisy  Channel  Algorithm.  The  important 
argument  that  the  S2^(x)  state  leads  to  more  efficient  error  detection 
than  the  S9n(x)  state  is  supported  by  the  fact  that  t2Q(x)  >  t01 (x) 


20 

unless  P12  »  PQ2. 


'21 


t22(x)  exhibits  some  interesting  behavior.  Figure  8  is  a  graph 
of  T22(x),  at  optimal  policy,  for  the  case  Pn0  =  P, 0  *  .1.  Note  that 


02  12 


the  abscissa  is  logarithmic  to  the  base  ,  and  therefore  increasing  x 
corresponds  leftward  motion  on  the  axis.  t22(x)  is  strictly  defined 
only  at  discrete  values  of  x.  The  function  reaches  a  local  maximum 
a  x  =  x£,  and  this  appears  to  be  a  general  property  of  t2^(x)  at  the 


optimum  policy  for  a  given  P. 


,U3 


As  x  increases  from  T(^)  toward  T,  t22(x)  decreases  at  first. 
The  time  to  process  a  S^Cx)  interval  does  not  always  increase  with 
the  length  of  the  interval.  This  can  be  explained  by  the  fact  that  the 
Noisy  Channel  Algorithm  policy  for  S^2  intervals  smaller  than  the 
cutoff  length  is  highly  inefficient.  Conflict  intervals  containing  _>  2 
messages  are  being  transmitted  in  their  entirety.  As  we  move  away  from 
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Capacity  for  P  -  0 


:  P  p 

02’  12 

NCTOFF* 

* 

T 

* 

X 

0 

10 

1.265 

.487 

10-6 

10 

1.265 

.487 

10“  3 

10 

1.270 

.485 

CN 

1 

O 

H 

9 

1.305 

.471 

2xl0“2 

6 

1.330 

.458 

5xlO~2 

5 

1.430 

.424 

7. 5x10“  2 

4 

1.460 

.402 

io_1 

3 

1.425 

.383 

TABLE  2 


(x)  at  Optimal 


NCTOFF  =  NCTOFF*= 


the  cutoff  point  there  is  a  chance  that  the  transmission  at  a 

state  will  result  in  a  message  being  successfully  transmitted  and  a 

S^(-)  interval  with  exactly  one  message  remaining.  This  is  then 

processed  quickly  and  efficiently.  Eventually,  as  we  move  further  away 

from  xc>  the  probability  of  an  interval  containing  more  than  two 

messages  grows,  and  this  results  in  t^C*)  a8ain  increasing  with  x. 

Thus  far  only  the  case  where  P  “  P^2  ^as  ^een  analyzed. 

Table  3  lists  the  optimal  algorithm  parameters  and  capacity  for  a 

number  of  asymmetric  cases,  where  and  P^2  are  unequal.  The  two 

corresponding  symmetric  cases  are  also  given  for  each  related  pair  of 

asymmetric  probabilities.  An  upper  limit  of  ten  was  imposed  on 
*  * 

NCTOFF  ,  and  T  was  optimized  with  a  grid  size  of  0.05. 

*  * 

For  notational  convenience  we  explicitly  indicate  X  *  T  ,  and 
* 

NCTOFF  as  functions  of  and  p-j.2  * 

* 


a  (p02.  p12) 


T  (P02,P12) 


NCTOFF  (P02»P12) 


A  few  behavioral  tendencies  of  the  assymmetric  error  probability 
cases  can  be  extrapolated  from  the  numbers  in  table  3.  First  of  all 
we  note, 


A  (P02,P12)  ~  *  (P12,P02) 


That  is,  Pq2  and  P^2  seem  to  have  roughly  the  same  influence  in  down¬ 
grading  channel  capacity  (capacity  seems  to  fall  off  slightly  faster 


with  P12) . 
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* 

* 

* 

P02 

P12 

NCT0FF 

T 

A 

.001 

.001 

10 

1.25 

.485 

.1 

.1 

3 

1.45 

.383 

.001 

.1 

6 

1.25 

.420 

.1 

.001 

4 

1.50 

.429 

0.0 

0.0 

10 

1.25 

.487 

0.1 

0.1 

3 

1.45 

.383 

0.0 

0.1 

6 

1.25 

.421 

0.1 

0.0 

4 

1.50 

.430 

0.0 

0.0 

10 

1.25 

.487 

o 

to 

0.2 

2 

1.55 

.323 

0.0 

0.2 

4 

1.30 

.370 

0.2 

0.0 

3 

1.65 

.391 

0.1 

0.1 

3 

1.45 

.383 

0.2 

0.2 

2 

1.55 

.323 

0.1 

0.2 

2 

1.30 

.344 

0.2 

0.1 

2 

1.45 

.353 

0.2 

0.2 

2 

1.55 

.323 

0.4 

0.4 

0 

1.30 

.219 

0.2 

0.4 

1 

1.30 

.260 

0.4 

0.2 

1 

1.70 

.271 

0.0 

0.0 

10 

1.25 

.487 

0.5 

0.5 

0 

1.60 

.151 

0.0 

0.5 

2 

1.25 

.247 

0.5 

0.0 

2 

2.65 

.259 

TABLE  3 
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In  the  regions  considered  here,  the  average  of  the  two  related 


symmetric  capacities  seems  to  be  a  fairly  good  approximation  to  an 
asymmetric  capacity,  or  "i 

^  ^02*^12^  +  ^ 


02*  02' 


12  12 


The  asymmetric  capacity  is  in  general  a  little  smaller  than  this 
average.  The  approximation  is  better  for  PQ2>  P  than  for  PQ2  <  P^. 
In  the  table  3  examples  it  is  also  evident  that. 


for  P 


max 


max  [P02,P12] 


P  .  =  min  [P„,P, „  ]  , 

mm  02  12  J  * 


NCTOFF* (Pmav »Pmav)  1NCTOFF*(P02,P12)  <  NCTOFF* (Pmin.Pmin>  •  No 

* 


max  max  — 

statement  can  be  made  for  T  . 

In  cases  where  P02  >  P12>  the  optimal  policy  calls  for  a 
and  a  smaller  (or  equal)  size  NCTOFF  than  would  be  optimum  for 
system  created  by  interchanging  P^  and  P^.  In  other 
words. 


equivalent 

larger  T 
the 


NCTOFF*(P02,P12)  <  NCTOFF* (P12,P02) 

T  (P02,P12)  -  T  (P12’P02)* 

This  is  understandable.  As  Pq2  becomes  relatively  large  in  comparison 
to  P^2*  becomes  desirable  to  increase  the  expected  number  of  messages 
in  a  Sq  transmission  interval.  A  P^2  error  can  only  occur  if  the 
transmission  interval  contains  no  messages  and  increasing  T  decreases 
the  probability  of  this  happening. 


5.4  Suggestions  for  Further  Research 


There  are  several  areas  of  further  study  suggested  by  this 
thesis  that  could  prove  interesting.  The  behavior  of  the  intermediate 
functions  derived  in  the  Noisy  Channel  Algorithm  could  be  analyzed  more 
rigorously.  The  system  behavior  for  asymmetric  error  probabilities  could 
also  be  characterized  in  depth.  There  are  also  more  general  issues 
concerned  with  introducing  errors  into  a  contention  resolution  scheme. 
How  does  the  introduction  of  a  wider  class  of  errors  effect  the 
functioning  of  the  Noisy  Channel  Algorithm?  How  can  a  transmission 
interval  algorithm  be  refined  to  more  effectively  process  different 
errors? 


i 


APPENDIX  A 


State  Transition  for  the  Noisy  Channel  Algorithm 

In  this  appendix  the  state  transitions  that  define  the  Noisy 
Channel  Algorithm  are  presented,  "q"  is  used  to  represent  the  interval 
that  is  on  the  top  of  the  stack.  When  the  next  state  is  indicated  as 
S(q),  this  implies  an  interval  was  taken  off  the  stack. 

Although  the  policy  states  (S2.  ^20^  an<*  ^1  P°^cy 

states  (S^,  S^q)  do  not  map  directly  to  the  Gallager  Algorithm  S2  and 
policies  the  lag  and  interval  length  are  updated  in  the  same  manner,  as 
implied  by  the  state  transition.  A  transmission  resulting  in  an 
idle  channel  being  detected  is  updated  as  follows: 


d  *  d,,  +  1  -  £  ,  , 

new  old  old 


£  =  £ 
new  o 


A  "transition"  from  a  S2  policy  state  to  S^(x)  results  in: 


d  =  d 
new  old 


new 


2£ 


old 


The  lag  and  interval  length  following  a  S2  state  are  updated 

It  it 

as  they  would  be  for  a  S2  state.  An  S^(x)  transition  to  S2(x)  implies, 

d  -  d  +  1 
new  old 


old 


new 
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Note  that  when  an  interval  is  taken  off  the  stack 


new 


1  +  d  i  .. 

old  old 


£  ■  q 

new 

The  state  classes  can  be  summarized  as  follows: 


Regular  States 


Policy 


0 

1 

2 


S 

S 

s 


0 

1 

2 


Error  States 


S 

S 

S 


10 

21 

20 


Policv 


S 

S 

s 


1 

2 

2 


Special  Processing  States 

* 

S1 

* 


Policy 

* 
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Regular  States 


Sent 

Detected 

Next  State 

V 

0 

0 

so 

0 

1 

S20(T> 

1 

1 

so 

1 

>  2 

S21(T> 

1  2 

I  2 

S2(T) 

Sx(x) 

:  0 

0 

won*t  occur 

0 

2 

•i  it 

1 

1 

so 

1 

1  2 

S21(x) 

>_  2 

>  2 

S2(x) 

S2(x) 

:  x  >  x  0 

c 

0 

S2^2)2 

0 

>  2 

s20(f>i  add 

(-) 

'2'2 

to 

stack 

1 

1 

S1^2 

1 

1 

S21(f)l  add 

(— ) 
k2'2 

to 

stack 

>  2 

* 

1  2 

S2(f)1  add 

(— ) 

2'  2 

to 

stack 

S  (x)  :  x  <  x  =+  S* (x) 
l  —  c  1 


Next  State 


S20(2^2 

S20<2)1  3dd  $2 

won't  occur 

If  It 

M  It 


S21<§)2 

s2o(!>i  add  <!>: 
S10(2)2 

S2i(2)i  add  (2^: 
won't  occur 

S1Q(q)  if  0  in  ( 
S1 (q)  if  >1  in 
Sq  if  stack 

S20(x) 


to  stack 


to  stack 

to  stack 


q 

empty 


won't  occur 


Special  Processing  States 


Sent 

Detected 

Next  State 

0 

0 

S20(q) 

if  0  in  q 

s2i<q> 

if  1  in  q 

S2(q) 

if  >  2  in  q 

so 

if  stack  empty 

0 

>2 

S*(x) 

1 

1 

S10(q) 

if  0  in  q 

Sx(q) 

if  >1  in  q 

so 

if  stack  empty 

1 

>2 

S*(x) 

>2 

>2 

S*(x) 

0 

0 

S20(2)2 

0 

>2 

S20('2)l 

add  (21)  to  stack 
2  2 

1 

1 

CM 

✓—s 

w 

O 

i— ♦ 

C/5 

if  0  in  (|)2 

S!<!>2 

if  >  1  in  (|)2 

1 

>2 

S21(I>1 

add  (y)2  to  stack 

S„(4),  add  (4)„  to  stack 


2 


>2 


APPENDIX  B 


Linear  Approximation  to  t22(x)  anc*  R(x) 

For  both  t 22 (x)  and  R(x)  we  need  to  make  an  approximation  to  the 
function  at  x  -  0.  This  then  yields  the  initial  points  necessary  to 
recursively  compute  t22(T)  and  R(T). 

In  linearizing  both  functions  Taylor  series  approximations  are  used 
on  the  Poisson  probability  distributions.  For  x  =  0,  the  following 
approximations  are  generated: 


1)  Pr (0  in  ||  >  2  in  x) 


-x/2 -x/2  x  -x/2) 
e  (1  -  e  -  —  e 

,  -x  -x 
1  -  e  -  xe 


5x_ 

48  „  1 


2)  Pr (1  in  ||  _>  2  in  x) 


X 

24 


3)  Pr(>^  2  in  >_  2  in  x) 


4)  Pr (1  in  ||  >1  in  |) 


5)  Pr  (  >  2  in  ||  >_  1  in  |) 
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this  leads  to, 


a  =  0 

b  =  1/6 

When  we  think  of  R(x)'s  interpretation  it  makes  sense  that 

lim  R(x)  =  0.  Nothing  can  be  returned  to  the  waiting  interval  from 
x-K) 

an  interval  of  zero  length. 

In  the  computer  programs  that  have  been  used  to  calculate  the 

capacity  the  linear  approximations  were  made  at  a  point  x  <_  T  •  2  ^ 

-9 

If  linear  approximations  at  x  =  2  are  compared  with  a  recursive 
-9-9 

calculation  of  R(2  )  or  t ^ (2  ),  based  on  an  approximation  at  x  =  2 

the  results  are  within  .01%  of  each  other.  This  accuracy  is  a  result 
of  the  probability  distributions  approaching  their  limiting  values. 


APPENDIX  C 


Capacity  Calculation  Computer  Programs 

Two  FORTRAN  computer  programs,  that  were  used  in  calculating  the 
Noisy  Channel  Algorithm  capacity  at  different  error  levels,  are  included 
in  this  appendix. 

The  first  program  determines  the  optimum  policy,  for  a  given 

* 

Pq2  and  over  a  specified  range  of  T  and  NCTOFF  values.  The 

second  program  performs  the  same  calculation,  but  it  sets  P^9  «  P^  and 

determines  the  optimum  policy  for  a  number  of  error  probability  values. 

Thus,  the  second  program  is  helpful  in  analyzing  how  the  algorithm 

behaves  in  relation  to  the  probability  of  error.  The  flexibility  of 

considering  asymmetric  cases,  however,  is  lost. 

In  these  programs,  the  linear  approximations  to  R(x)  and  t^^(x) 

takes  place  at  x  <  2  If  x  >  2  ^  then  the  t„„  recursion  formula 

—  c  22 

of  equation  4.9  is  used  to  find  t^Cx^,  and  the  formula  of  equation  4.10 
is  then  used  to  determine  t^^lT).  The  same  R(x)  recursion  formula  of 
equation  4.12  is  used  for  any  value  of  x.  This  procedure  of  starting 
recursion  at  x  <  2  ^  insures  good  approximations  to  the  two  functions 
(as  discussed  in  Appendix  B) .  It  also  insures  that  the  computational 
accuracy  of  the  program  is  relatively  independent  of  the  algorithm 
parameters  being  investigated. 

Double  precision  variables  are  used  in  both  programs.  The  programs 
begin  to  become  significantly  inaccurate  as  T  approaches  zero  or  as 
NCTOFF  approaches  30.  In  these  cases,  the  interval  lengths  being 


considered  are  extremely  small  and  some  of  the  distribution  function 


calculations  suffer  from  important  round-off  effects.  It  has  been 
mentioned  that,  as  far  as  algorithm  efficiency  is  concerned,  NCTOFF  =  10 
is  a  good  approximation  to  any  larger  NCTOFF  value. 

The  maximum  stable  throughput  achievable  for  a  given  system  appears 
to  be  a  convex  function  of  NCTOFF.  For  a  given  NCTOFF,  the  maximum 
throughput  appears  to  be  a  convex  function  of  T.  These  two  properties 
combine  to  make  it  relatively  simple  to  find  the  optimal  policy 
corresponding  to  a  set  of  error  probabilities. 

In  running  these  programs  the  t2o’t21,t22  °r  R  functions  can  be 
studied  by  strategically  inserting  a  few  WRITE  statements. 
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THE  NOISY  CHANNEL  ALGORITHM 
This  program  optimizes  system  throughput  over 
a  designated  range  of  T  and  NCTOFF  values. 

IMPLICIT  REAL*8  (A-H.O-Z) 

EXP(X)  =  DEXP(X) 

10  WRITE (6,) 'INPUT' 

Input  error  probabilities  and  boundaries  for 
T  and  NCTOFF  values  considered.  NCTOFF  is 
considered  at  integer  values,  and  T  is 
optimized  at  values  spaced  by  TGRID.  1=1 
indicates  that  maximum  throughput  at  all 
NCTOFF, T  pairs  is  to  be  outputted. 

Otherwise,  only  optimum  is  printed  out. 

INPUT, P02.P1 2,TL, TH ,TGRID, NCTOFL, NCTOFH, I 
If  TL=0  then  exit  program. 

IF  (TL  .EQ.  0.)  GOTO  111  0 
Compute  t20*  and  t21 *  . 

T20CUT= (1.0  +  P0  2 ) / ( 1 .0-P02-P02*P02) 

T21CUT=(1 .0+1 . 5  *  P 1 2+P1 2*P02*T20CUT ) /( 1 ,0-( PI 2*P1 2+P1 2) /2 .0) 
Compute  t22  linearization  coefficients. 

A=5 .0+P02"T20CUT/2 .0+2 .0"P1 2*T21 CUT 
B=( 1 9 . 0+P02*T20CUT+3 . 0«P 1 2«T2 1  CUT ) /I  8 . 0 
NUMBER  =  number  of  T  iterations  necessary. 

NUMBERr(TH-TL) /TGRID  +1.6 

0CAP,0T,  and  NCTOF  indicate  optimum  values. 

OCAP=0 .0 

Loop  over  T  values. 

DO  350  K=1, NUMBER 

Increment  T  by  TGRID  at  beginning  of  new  loop 
through. 

T=TL+(K-1) "TGRID 

Loop  over  NCTOFF  values. 

DO  350  NCT0FF=NCT0FL, NCTOFH 

L  is  the  number  of  divisions  beyond  cutoff 

point  necessary  to  insure  linearization  accuracy. 

L=0 

50  IF  (T"( ,5**( NCTOFF  +  L) )  .LE.  ( .5**10) )  GOTO  90 
L=L+1 
GOTO  50 

LIMIT  is  the  total  number  of  recursions  from 
the  linearization  point  to  T. 

90  LIMIT  =  NCTOFF  +  L 

Make  the  linear  approximations  to  t22  and  R. 

T22CUT  =  A  +  B»T»( .5**LIMIT) 

X  =  T»(.5”LIMIT) 

RETURN  =  X/6.0 
T22=T22CUT 
T21=T21CUT 
T20=T20CUT 

Begin  recursion  towards  T. 

DO  300  J=1, LIMIT 

Double  interval  length  each  time  through. 

Y  =  new  length;  X  =  old  length. 

Y=2.0»X 
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C  Compute  conditional  Poisson  probability  distributions. 

DEN0M1 =1 . 0-EXP (-Y)-Y*EXP(-Y) 

P0ISS1 =(EXP( -X)*( 1 . 0— ( 1 .0+X)*EXP(-X) ) )/DEN0M1 
P0ISS2=(X*EXP( -X ) *( 1 . 0-EXP ( -X ) ) )/DEN0M1 
P0ISS3=( 1-(X+1 ) *EXP{ -X) )/DEN0M1 
DEN0M2=X*EXP( — X ) * ( 1 .0-EXP( -X) ) 

P0ISS5=((X»EXP(-X))»(1 .0-(1 .0+X)»EXP(-X)))/DEN0M2 
P0ISS4=(X  *EXP( -X) )/( 1 .0  -  EXP(-X)) 

C  R  recursion  formula. 

RETURN  =  RETURN  •  (1.0  -  POISS4«POISS2)  +  X*P0ISS3 
If  we  are  below  cutoff  then  we  only  do  t22 
recursion  before  cycling  back.  Otherwise  go  to 
150  and  do  all  three  recursions. 

IF  ((J  -  L)  .GE.  1)  GOTO  150 

T22CUT  =  2.0+P0ISS1 *(P02*T20CUT+T22CUT)  +  P0ISS2»(P12»T21CUT) 

&*(1.0  +  P0ISS4 )  +  1.0  +  P0ISS5*T22) )  +  P0ISS3»(1.0  T22) 

T22  =  T22CUT 
GOTO  290 

150  T22  =  (P0ISS1  *(  1  .0+P02*T20+T22)  )+( P0ISS2*( 2.0+P1 2*T21* ( 1.0  +P0ISS14) 
&+POISS5»T22))+POISS3*(1.0  +  T22/ 

T21 si  .5  +  P02#T20+( 1.0  +P12)«T2l/2.0 
T20s1 .0  +  (1 .0+P02) *T20 
Update  X. 

290  X=2.0»X 
300  CONTINUE 

Compute  Poisson  probabilities  for  0,1, >=2  messages  in  T. 
PR0B0=EXP(-T) 

PR0B1 =T#PROB0 
PR0B2=1 .0-PR0B1 -PROB0 

Compute  the  maximum  achievable  stable  output  for  this  set  of 
T  and  NCTOFF  values. 

TOP  =  T-PR0B2* RETURN 

CAPCTY  =  T0P/( 1.0  +  (PROB0»P02»T20)  ♦  ( PRC® 1 *P1 2*T21 )  +  (PR0B2»T22>) 
IF  (I  .EQ.  1)  WRITE(6,)T, NCTOFF, CAPCTY 
Have  we  found  a  new  optimum  policy? 

If  not.  cycle  back  Immediately  for  new  NCTOFF, T  pair. 

If  we  do  have  a  new  maximum  throughput,  the  optimum 
parameters  and  throughput  are  updated  first. 

IF  (CAPCTY  ,LT.  OCAP)  GOTO  350 
OCAP  =  CAPCTY 
NCTOF  =  NCTOFF 
OT  s  T 
350  CONTINUE 

Output  maximum  throughput  and  optimal  policy  for  the 
range  of  algorithm  parameters  that  have  been 
considered.  Loop  back  for  new  input. 

WRITE ( 6 ,) OT , NCTOF , OCA  P 
1)00  GOTO  10 
410  STOP 
END 


n  n  no  o  o  non 


C  THE  NOISY  CHANNEL  ALGORITHM 

C  This  program  is  virtually  identical  to  the  previous 

C  program  that  finds  the  optimal  policy,  for  a  given  P02 

C  and  PI 2,  over  a  range  of  NCTOFF  and  T  values.  In 

C  this  program,  however,  we  assume  P02=P12,  and  the  optimal 

C  policy  is  computed  over  a  range  of  error  probabilities. 

C  The  same  NCTOFF  and  T  limits  are  used  in  examining  the 
C  policy  at  different  error  probabilities.  Thus  this 

C  program  just  involves  adding  an  outer  DO  loop  to  the 

C  the  previous  program.  This  loop  just  cycles  through 

C  different  error  probabilities.  We  therefore  can  then 

C  observe  how  the  optimal  policy  and  capacity  vary  with 

C  channel  characteristics. 

IMPLICIT  REAL*8  (A-H.O-Z) 

EXP ( X )  =  DEXP(X) 

10  WRITE(6,) 'INPUT1 

PL  and  PH  are  the  lower  and  upper  limits  respectively 
on  the  error  probabilities  that  are  considered.  PGRID 
is  the  spacing  of  the  probabilities  . 

INPUT , PL , PH , PGRID , TL , TH , TGRID , NCTOFL , NCTOFH 
IF  (TL  .EQ.  0.)  GOTO  410 

NUMB  is  the  number  of  iterations  with  different 
probabilities  to  be  done. 

NUMB= ( PH- PL) /PGRID  +  1.6 

Now  establish  the  external  Do  loop  that  cycles  the 
program  with  different  error  parameters. 

DO  390  M= 1 , NUMB 

Update  the  error  probabilities  at  the  beginning  of 
each  loop  through. 

P02=PL+PGRID*(M-1 ) 

PI 2=P02 

T20CUT= (1.0  +  P02)/(1 ,0-P02-P02*P02) 

T21CUT=( 1 .0+1 .5*P12+P12*P02*T20CUT)/( 1 .0-(P1 2*P1 2+P12)/2.0) 
A=5.0+P02»T20CUT/2.0+2.0»P12»T21CUT 
B= ( 1 9 . 0+P0  2»T2  0CUT+3 • 0  *  P 1 2»T2 1 CUT ) / 1 8 . 0 
NUMBERS (TH-TL) /TGRID  +  1.6 
OCAPsO.O 

DO  350  Ksl .NUMBER 
TsTL+(K-1 ) *TGRID 
DO  350  NCTOFFsNCTOFL, NCTOFH 
LsO 

50  IF  (T*(.5**(NCT0FF  +  L) )  .LE.  (.5”10))  GOTO  90 
LsL+1 
GOTO  50 

90  LIMIT  s  NCTOFF  +  L 

T22CUT  s  A  +  B»T»( .5»*LIMIT) 

X  s  T«(.5#*LIMIT) 

RETURN  s  X/6.0 
T22sT22CUT 
T21 sT21CUT 
T20sT20CUT 
DO  300  Jsl, LIMIT 
Ys2.0»X 

DENOMIsl .  0-EXP (-Y)-Y»EXP(-Y) 
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P0ISS1 =(EXP(-X)*( 1 .0— ( 1 .0+X)#EXP(-X) ) )/DENOM1 
POISS2 = ( X»EXP ( -X ) • ( 1 . 0-EXP ( -X ) ) ) /DENOM 1 
POISS3 = ( 1 -( X+1 ) «EXP ( -X ) ) /DENOM 1 
DENOM2 =X*EXP ( -X )•( 1 . 0-EXP ( -X ) ) 

POISS5  =  ((X»EXP(-X))»(1 .0-0  .  0+X ) *EXP ( -X ) ) ) /DEN0M2 
P0ISS4=(X  *EXP(-X) )/( 1 .0  -  EXP(-X)) 

RETURN  =  RETURN  «  (1.0  -  P0ISS4*P0ISS2)  +  X»P0ISS3 
IF  ((J  -  L)  .GE.  1)  GOTO  150 

T22CUT  =  2.0+P0ISS1 *(P02*T20CUT+T22CUT)  ♦  P0ISS2»(P1 2*T21CUT) 

*•(1.0  +  P0ISS4 )  +  1.0  +  P0ISS5*T22) )  +  P0ISS3«(1.0  +  T22) 

T22  =  T22CUT 
GOTO  290 

150  T22=(P0ISS1 •( 1 .0+P02*T20+T22) )+(POISS2*(2.0+P12*T21B(  1 .0  +P0ISS4) 
&+P0ISS5*T22) )  +  POISS3«(1.0  +  T22) 

T21=1.5  +  P02#T20+( 1.0  +P12)»T21/2.0 
T20=1 .0  +  ( 1 .0+P02) *T20 
290  X=2.0»X 
300  CONTINUE 

PROBOsEXP(-T) 

PR0B1 =T*PROBO 
PR0B2=1 .0-PR0B1 -PROBO 
TOP  s  T-PR0B2«RETURN 

CAPCTY  r  TOP/O.O  +  ( PROB0#PO2*T2O)  ♦  ( PRCB1  «P1 2»T21 )  +  (PR0B2*T22)) 
IF  (CAPCTY  .LT.  OCAP)  GOTO  350 
OCAP  =  CAPCTY 
NCTOF  =  NCTOFF 
OT  =  T 
350  CONTINUE 

C  The  maximum  throughput  and  optimal  policy  are  printed 

C  out  at  each  error  probability  considered. 

WRITE( 6 , ) P02 . OT, NCTOF, OCAP 
390  CONTINUE 
400  GOTO  10 
410  STOP 
END 


References 


Capetanakis,  J. ,  The  Multiple  Access  Broadcast  Channel: 

Protocol  and  Capacity  Considerations,  Ph.D.  Thesis,  M.I.T.  ,1977. 


Drake,  Alvin  W, ,  Fundamentals  of  Applied  Probability  Theory. 
McGraw-Hill,  Inc.,  1967. 

Mosely,  J.,  An  Efficient  Contention  Resolution  Algorithm  for 
Multiple  Access  Channels,  M.S.  Thesis,  M.I.T. ,  1979. 


Distribution  List 


Defense  Documentation  Center  12  Copies 

Cameron  Station 
Alexandria,  Virginia  22314 

Assistant  Chief  for  Technology  1  Copy 

Office  of  Naval  Research,  Code  200 
Arlington,  Virginia  22217 

Office  of  Naval  Research  2  Copies 

Information  Systems  Program 
Code  437 

Arlington,  Virginia  22217 


Office  of  Naval  Research 
Branch  Office,  Boston 
495  Summer  Street 
Boston,  Massachusetts  02210 

Office  of  Naval  Research 
Branch  Office,  Chicago 
536  South  Clark  Street 
Chicago,  Illinois  60605 

Office  of  Naval  Research 
Branch  Office,  Pasadena 
1030  East  Greet  Street 
Pasadena,  California  91106 

New  York  Area  Office  (ONR) 

715  Broadway  -  5th  Floor 
New  York,  New  York  10003 

Naval  Research  Laboratory 

Technical  Information  Division,  Code  2627 

Washington,  D.C.  20375 

Dr.  A.  L.  Slafkosky 
Scientific  Advisor 

Commandant  of  the  Marine  Corps  (Code  RD-1) 
Washington,  D.C.  20380 


1  Copy 


1  Copy 


1  Copy 


1  Copy 


6  Copies 


1  Copy 


Office  of  Naval  Research 
Code  455 

Arlington,  Virginia  22217 


1  Copy 


Office  of  Naval  Research  1  Copy 

Code  458 

Arlington,  Virginia  22217 

Naval  Electronics  Laboratory  Center  1  Copy 

Advanced  Software  Technology  Division 
Code  5200 

San  Diego,  California  92152 

Mr.  E.  H.  Gleissner  1  Copy 

Naval  Ship  Research  &  Development  Center 
Computation  and  Mathematics  Department 
Bethesda,  Maryland  20084 

Captain  Grace  M.  Hopper  1  Copy 

NAICOM/MIS  Planning  Branch  (OP-916D) 

Office  of  Chief  of  Naval  Operations 
Washington,  D.C.  20350 


Advanced  Research  Projects  Agency  1  Copy 

Information  Processing  Techniques 
1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

Dr.  Stuart  L.  Brodsky  1  Copy 

Office  of  Naval  Research 
Code  432 

Arlington,  Virginia  22217 

Captain  Richard  L.  Martin,  USN 
Commanding  Officer 
USS  Francis  Marion  (LPA-249) 

FPO  New  York  09501 


1  Copy 


